Kinks: Fingerprints of strong electronic correlations 
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Abstract. The textbook knowledge of solid state physics is that the electronic specific heat 
shows a linear temperature dependence with the leading corrections being a cubic term due to 
D ' phonons and a cubic-logarithmic term due to the interaction of electrons with bosons. We have 

»' | , shown that this longstanding conception needs to be supplemented since the generic behavior 

of the low-temperature electronic specific heat includes a kink if the electrons are sufficiently 
strongly correlated. 



1. Introduction 

Landau's Fermi liquid theory pQ can be considered as the "standard theory" of solid state physics. 
It predicts that the electronic properties of a (normal) metal follow, in the presence of Coulomb 
correlations, a renormalized (quasiparticle) version of those for non-interacting electrons. Among 
■ the most fundamental physical properties is the heat capacity which, at low temperatures, is 

predominated by the electronic degrees of freedom. For a normal metal, Landau's Fermi liquid 
theory [T] predicts a linear increase of the specific heat capacity with temperature cy = jflT 
and a cubic term, more precisely a term ~ T 3 log(l/T), as the leading correction [2J [3]. For 
free electrons, the prefactor 70 is directly proportional to the electronic density of states at the 
Fermi level since, due to the Pauli principle, only these electrons contribute. In the case of a 
Fermi liquid, we merely need to introduce a quasiparticle renormalization factor < Zfl < 1 to 
account for the Coulomb interaction, which enhances the specific heat of a correlated metal by 
Ifl = 7o/-^fl- However, this electronic contribution prevails only at temperatures much lower 
than the Debye temperature. Otherwise the cubic phononic term, which has a much higher 
«-j ' prefactor because of its bosonic nature, dominates. 

As we have shown in a recent letter [3], this longstanding conception needs to be 
supplemented: For strongly correlated electrons, a clear kink in the temperature dependence 
of the specific heat appears, marking the abrupt change from one linear behavior to a second 
one with a reduced slope at higher temperatures. This can be shown numerically solving the 
Hubbard model [5] within dynamical mean field theory [BJ, using exact diagonalization as an 
impurity solver. Quantitatively the same results can also been obtained analytically, following 
the derivation by Abrikosov et al. [2] for the specific heat, taking the recently observed kinks in 
the self energy as a starting point [8]. 

The experimental confirmation of these theoretical results is somewhat problematic since 
typically, at the kink temperature, the phononic contribution to the specific heat is already 
dominating. An exception in this respect is L1V2O4 [9l [10], the first heavy Fermion system 



with (i-electrons, since the kink temperature is here particularly low. And, indeed, the precise 
inspection of the experimental data shows a kink in the correct temperature range [4]. In this 
Proceeding, we will review the previous theoretical results and include additional data. 



2. Kinks in the self energy 

In the presence of strong electronic correlations, the k-integrated spectral function A{uj) shows 
a typical three peak structure with a lower and an upper Hubbard band and a quasiparticle 
peak in between. Using the DMFT formula 

E(w) =u + n-l/G(u) - A(G(u)) (1) 

which relates Green function and self energy at frequency u>, one can show [8] that such a strongly 
correlated three peak structure necessarily entails a kink. Note that the hybridization function 
A(G(o;)) in Eq. (JTJ) is for a simple semi-circular density of states just A(G(ui)) = ^-G(oj) (D: 
half the bandwidth). For other lattices, A(G(u)) depends on the moments of the non-interacting 
density of states with a similar term proportional G(u) dominating. 

The kink argument is now as follows [8]: The first three terms of the r.h.s of Eq. (pQ) yield an 
almost linear frequency dependence for the real part of the self energy over the entire region of 
the central quasiparticle peak. The derivative hence gives a quasiparticle renormalization factor 

Z QP = [1 - dReEM/c^Uo]- 1 (2) 

The last (hybridization) term however yields an additional contribution to ReE(cj) which 
basically is proportional to ReG(w). This real part can be directly obtained from the 
spectral function of the central quasiparticle peak A(u>) = — ImG(u) via a Kramers-Kronig 
transformation. Inside the overall width of the quasiparticle peak, there is a turning point in 
A(u). Hence, ReG(u;) has a linear frequency dependence up to an maximum at the turning 
point A(u). After this maximum changes in ReG(w) are minor. Altogether this means that 
at low frequencies we have to add the slope obtained from the A(G(uj)) in Eq. ([2]) yielding an 
altogether smaller Fermi liquid renormalization factor Zfl at low energies. In between the two 
regimes Zfl and Zqp, there is a kink at a frequency u*. In summary we hence have: 

f (1 - -J—)oJ UJ <UJ* 
ReE(w) = \ H i\ L , * (3) 

with the constant —b = (^^7 — z^) 1 ^* providing for a continuous function. 

Before we turn to the specific heat in the next section, let us note that the kink in the self 
energy directly leads to a kink in the energy-momentum dispersion relation of the correlated 
electrons. Let us also remark here, that the result of generic kinks in the self energy of strongly 
correlated systems should be robust beyond DMFT, as similar effects are to be expected in 
cluster [11] and diagrammatic extensions [12] of DMFT. 



3. Analytical formula for the specific heat 

Based on Eq. Q, we have developed an analytical theory for the specific heat on the basis of a 
formula by Abrikosov, Gor'kov and Dzyaloshinski (AGD) for the entropy of a fermionic system 
at low temperatures [2]: 

1 f°° f°° e y 

S(T) = —J ^deN(e) J ^ dyy - [logG R \e,yT) - \ogG A \e,yT)l (4) 

This AGD formula is (through a low-temperature diagrammatic expansion) based on the low 
frequency behavior of the self energy E(u;) at zero temperature [or the corresponding retarded 



and advanced Green function in Eq. ([4"D] so that we can directly apply it with from 
Eq. ([3]). Let us note that iV(e) is the density of states and the frequency integral has been 
rewritten through a dimensionless variable y = uj/T (kg = 1). The AGD formula can 
be easily differentiated w.r.t. T, so that the specific heat cy(T) = T dS ^ can be computed 
directly. With some algebra (see [3] for more details) and substituting substituting A(e, yT) by 
5-functions, we obtain the final result 
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Knowing the two renormalization factors Zpp and -Zqp and the kink frequency uj* as well 
as the non-interacting density of states N(e), Eq. ([5]) allows us to calculate the specific heat 
through a simple integral. If vice versa, the specific heat is known, e.g., experimentally, we can 
employ Eq. ([5]) as a fit formula for obtaining Zpp, Zqp and ui* . 



4. Experimental validation 

The latter (fit formula) approach had been taken for the experimental validation of the specific 
heat kinks in L1V2O4 since LDA+DMFT calculations [13J for L1V2O4 [T3] could not provide for 
a fine enough resolution of the low frequency dependence of the self energy, even when using 
the projective quantum Monte Carlo method [TO] which is most appropriate for this purpose. 
As Fig. 3 of [12] shows the three fitting parameters provide for an excellent agreement of Eq. 
© and the specific heat of [91 [TO], including a clearly visible kink. Let us note that we did 
not subtract a phonon contribution for the specific heat since (i) it does not alter the existence 
or non-existence of a kink, (ii) it is still an order of magnitude smaller than the electronic 
contribution at the kink temperature and (iii) there is a considerable arbitrariness depending 
on the reference system chosen. 



5. Numerical results 

Fig. [1] shows results of our final equation ([5]) applied to the Hubbard model with a semi-circular 
density of states with half bandwidth D = 1. As an input to Eq. ([5]), Zpp, Zqp and ui* have been 
fitted to the numerical renormalization data for the DMFT self energy from Ref. [TO]. Besides, 
we have also calculated the DMFT specific heat directly, using exact diagonalization (ED) with 
7 levels in the impurity bath. Both results agree very well, and the minor deviations can be 
attributed to the numerical inaccuracy of the exact diagonalization study for which cy had to 
be obtained via a numerical differentiation of the total energy. 

In the inset, we show the specific heat over a slightly extended temperature range. As the 
quasiparticle peak eventually ends, the specific heat has a maximum after the kink temperature. 
Note that this maximum is not contained in our formula © since it would require an extension 
of the self energy description Eq. (|3|) to higher energies where ReS(u;) has a maximum at the 
frequencies between Hubbard and quasiparticle bands and ImS(w) becomes important. 



6. Conclusion 

Landau's Fermi liquid theory predicts a linear temperature dependence of the electronic specific 
heat capacity with temperature. We have shown this longstanding conception needs to be 
supplemented: For strongly correlated electrons, a clear kink in the temperature dependence 
appears, marking the abrupt change from one linear behavior to a second one with a reduced 
slope at higher temperatures. Recent experiments on L1V2O4, an ideal material for studying 
the electronic specific heat, confirm our theory. A consequence of our findings is that materials 
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Figure 1. Specific heat as a function of temperature for various values of the Coulomb 
interaction U obtained via the AGD based equation ([5]) [solid lines] and numerically by exact 
diagonalization of the DMFT equations (crosses). Inset: behavior over an extended temperature 
range which shows a maximum after the kink temperature marking the end of the quasiparticle 
band. 

with correlated electrons are more resistive to cooling at low temperatures (where cooling is 
particularly difficult) than expected from the behavior at intermediate temperatures. 
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